COMBINATORIAL MINORS OF MATRIX FUNCTIONS 
AND THEIR APPLICATIONS 



■ VLADIMIR SHEVELEV 

^~>, Abstract. For every matrix function we introduce its "combinatorial 

minors" . As important examples of their applications, we give a recur- 
sive algorithms for solution of several general enumerative problems of 
permutations with restricted positions. 



1. Introduction 



o 
u 

1^ . Let iS'-"^ be symmetric group of permutations of numbers {1, ...,n}. Let 

A = {O'ij} be square matrix of order n. Recall that permanent of A is defined 
by formula (^) 



>. (1-1) perA= ^ JJa,,,(,). 

in 
cn 



i=l 

If y4 is a (0, 1) matrix, then it defines a class B = B{A) of permutations with 
1^ - restricted positions, such that the positions of its zeros are prohibited. Such 

O ■ class could be equivalently defined by a simple inequality: a permutation 

TT e S if and only if for its incidence matrix P we have P < A. One of the 
most important application of per A consists of the equality \B\ = per A. 
^ I Thus per A enumerates permutations with restricted positions of the class 

Let 7(7r) be number of independent cycles of tt, including cycles of length 
1. Then the difference (i(7r) = n — is called decrement of n ([3]). Per- 
mutation 71 is called even (odd) if di^Tc) is even (odd). Note that determinant 
of matrix A one can define by the formula 

n n 

(1.2) detA= ^ n^*'^»~ X] n^'-^w- 

even se5(") «=1 odd sG5(") «=1 

Since, evidently, we also have 

n n 

(1.3) perA= ^ J] 

even sG5(") «=1 odd sG5(") «=1 

then the numbers of even and odd permutations of class B{A) are given by 
vector 

(1.4) {^{perA + detA), ^{perA-detA)). 
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Note that, in the contrast to permanent, there exist methods of very fast 
calculation of detA. Therefore, the enumerative information given by (II ■4p 
one can obtain approximately for the same time as the number \B\ given by 

Let m > 3 and < A; < m be given integers. We say that a permutation 
TT belongs to class k modulo m {n G S'^lj, if d{TT) = k (mod m). Let now A 
be (0, 1) square matrix of order n. The first problem under our consideration 
is the enumeration of permutations tt G i3 of class k modulo m. It is clear 
that this problem is a natural generalization of problem of enumeration of 
even and odd permutation with restricted positions which is solved by (1.4). 
In order to solve this more general problem, put u = and introduce a 
new matrix function which we call u-permanent. 

Definition 1. Let A be a square matrix of order n. We call u-permanent 
of matrix A the following matrix function 

n n 

i = l * = 1 



i,s{i)- 



j(n) i = l ^^^{n) i=l 



m — 1 ,m 

Note that, if m = 1, then per^iA = perA, and if m = 2, then per^A = 
detA. In case m > 3, every sum in (1.5) essentially differs from permanent. 
Therefore, the known methods of evaluation of permanent (pj, ch. 7) are 
not applicable. However, using so-called "combinatorial miners", below we 
find an expansion per^A over the first row of matrix A. This allows to re- 
duce a problem of order n to a few problems of order n — 1. 

The second problem under our consideration is another important prob- 
lem of enumeration of full cycles with restricted positions. In connection 
of this problem, we introduce another new matrix function which we call 
cyclic permanent. 

Definition 2. Let A be square matrix of order n. The number 

n 

(1.6) Cycl{A) = Y,X{(^^,s(i)^ 

s i=l 

where the summing is over all full cycles from we call a cyclic perma- 
nent of A. 

The third problem is a problems of enumeration of permutations with 
restricted positions with a restrictions on their cycle structure. 

DfiTintfi ^( thp, niimbpr nf all rvrlps nf nfirmiitatinn .s inrliidinp' rvrlfis nf 
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length 1. Recall ([I]) that the absolute value of Stirling number S{n,k) 
of the first kind equals to number of permutations s G S^"'^ with 7(5) = k 
(PP,[S])- ^ natural generalization of Stirling numbers of the first kind is the 
following matrix function. 

Definition 3. The matrix function 

n 

(1.7) S{A;n,k)= H'^^'^^W' 

where n is order of square matrix A, we call Stirling function of index k. 

Finally, recall (|S]) that a permutation s E 5'-"'^ with ki cycles of length 
1, ^2 cycles of length 2, and so on, is said to be of cycle structure k = 
{ki,k2, kn)- Denote iy{ki, ^2, ^n) the number of permutations of class 
k = {ki, k2, kn). Then polynomial 

(1.8) C(ti, t2, t„) = 5^ i^iki, k2, kn)t'lH'2'-t':: 

is called the cycle indicator of permutations of S^"^ . A natural generalization 
of the cycle indicator of permutations with restricted positions is 

(1.9) C{A; ti, t2, tn) = 5^ i^{A; k,,k2,..., kn)tlH',\ . .t^- 

where i^iA; ki,k2, kn) is the number of permutations of class B{A) with 
the cycle structure k. 

2. Observations in case m = 2 of cj-determinant 

The case of determinant (m = 2) is a unique case when it is easy to obtain 
a required enumerative information formally given by formulas (11. 4p . For 
the passage to a general case it is important for us to understand how one 
can obtain such information from the definition (11. 2p of determinant only. 
Essentially, the required information is contained in vector 

n n 

(2.1) detA={ ^ n^^'^W ' ~ 5Z n^^'''(^))' 

even s£S i=l odd sGcS i=l 

and it immediately disappears if to use an identity of the form 1 — 1 = 0. 
None of algorithms of fast calculation of determinant exists without such 
identity. On the other hand, in the Laplace algorithm of expansion of 
determinant over (the first) row such identity one can use only in the last 
step. Therefore, if not to do the last step, we can obtain the required 
enumerative information. 

Example 1. By the Laplace expansion, we have 

detll 1 1 =-(l-l) + (l-l) = -1 + 1 + 1- 1 = 2-2 
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and, if not to do the useless (with the enumerative point of view) last step, 
then we have 

_/0 1 1\ 
det\l 1 1 = (2,-2). 
\l 1 l) 

This means that there are two even and two odd permutations with the 
prohibited position (1,1). 

Unfortunately, most hkely, it is impossible to generalize for cu-permanent 
the Laplace expansion of determinant over the first row in its classic form 

n 

(2.2) detA = ^{-ly-^aijMij, 

i=i 

where Mij is minor of element a^, i.e., determinant of the complementary 
to aij submatrix Aij. 

Therefore, let us introduce a more suitable for our aims notion of a "com- 
binatorial minor" of element a^j. Let the complementary to aij submatrix 
Aij have the following n — 1 columns 

(2.3) Ci, C2, Cj_i, Cj+i, c„. 

The first j — 1 of these columns we change in the following cyclic order: 
C2, C3, Cj-i, c\. Then we obtain a new matrix Aij with the columns 

(2.4) 

Determinant of matrix Aij we call combinatorial minor {CM)ij of element 
aij . It is easy to see that 

(2.5) {CM)ii = Ma; (CM),, = {-ly-^Mij, j=2,...,n. 
Therefore, e.g., expansion (12.21) one can rewrite in the form 

n 

(2.6) detA = ai,i(CM)n - ^aij(CM)y. 

i=2 

In general, let us give a definition of combinatorial minors for arbitrary 
matrix function X{A). 

Definition 4. Let X be matrix function defined on all square matrices of 
order n > 3. Let A = {aij} be a square matrix of order n and Aij be 
the complementary to aij submatrix with columns (\2.!^ . Denote Aij a new 
square matrix of order n — 1 with columns ( |^.^p . Then the number X^A^j) 
is called a combinatorial minor of aij. 

It appears that our observation (12. 6 p has a general character. So, in 
Sections 4 we give a generalization of the Laplace expansion of type (12. 6 p 
for per^A, Cycl{A), S{A;n,k) and Cn(ti, ^2, ^n)- 



(3.2) vr*(^) 



(3.4) 7i*{i) 



combinatorial minors 5 
3. Main lemma 

Lemma 1. Let n E S''"'^ with = 1 and a = (Tj(7r), j > 2, such that 

a{l)=n{2), a{2) = n{3), a{j - 2) = ir{j - 1), a{j-l) = n{l), 

(3.1) a(j) = vr(j) = 1, a(i + l) = 7r(j + l), a(n) = 7r(n). 

Let, further, n* G S^""^^^ defined by the formula 

a(z)-l, if l<t<j-l, 
a{i + l) - 1, if j <i<n-l. 

Then permutations vr and n* have the same number of cycles: 

(3.3) 7(vr)=7(7r*). 

Proof. From (I3ll])-(l32!) we find 

'7r(^ + l)-l, z/ 
_7r(l) - 1, if i= j -1. 
Consider a cycle of n containing number j. Let it has length I > 2, such 
that 

7r(j) = 1, 7r(l) = fci, n{ki) = k2, nik^^) = ki_2, 'K{h-2) = 3- 

Beginning with the equality 7r(l) = ki, by (13. 4p . this means that 

7r*(j - 1) = A;i - 1, 7r*(A;i - 1) = A;2 - 1, 

7r*(A:,_3 - 1) = - 1, vr*(A;,_2 - 1) = J - L 
Thus to cycle of vr containing number j of length / > 2 corresponds a cycle 
of length / — 1 of TT*. Quite analogously, we verify that to cycle of vr not 
containing number j of length / > 2 corresponds a cycle of the same length 
of vr*. E.g., to cycle of length / > 2 of the form 

7r(j-l) = A;i, Txikx) = k2, 7r(/c2) = k^, rr{ki_2) = k-i, T^ik-i) = j-l 

(beginning with the equality n{ki) = k2), corresponds the cycle of the same 
length 

7r*(A:i - 1) = A;2 - 1, 7r*(A;2 - 1) = fcg - 1, 
7i*iki.2 - 1) = h^i - 1, 7r*{ki^i -l)=j-2 

such that - 2) = n{j - 1) - 1 = A;i - 1. 

■ 

Note that the structure of Lemma [1] completely corresponds to the pro- 
cedure of creating the combinatorial minors. 
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4. Laplace expansions of type (2.6) of per^A, Cycl{A), S{A;n,k) 

AND C{ti,t2, ...,tn) 

1) peVi^jA. Consider all permutations vr with the condition 7r(j) = 1. Let j 
corresponds to j-th column of matrix A. Then the considered permutations 
correspond to diagonals of matrix A having the common position (1, j). If 
j = 1, then, removing the first row and column, we diminish on 1 the number 
of cycles of every such permutation, but also we diminish on 1 the number 
of elements of permutations. Therefore, the decrement of permutations 
does not change. If j > 2, consider continuation of these diagonals in the 
matrix of combinatorial miner Aij. Then, by Lemmad], the number of cycles 
of every its diagonal does not change and, consequently, the decrement is 
diminished by 1. This means that we have the following expansion of per^^A 
over the first row 

n 

(4.1) det^A = ai,i(CM)n +a;^aij(CM)i,-, 

i=2 

where {CM)ij, j > 1, are combinatorial minors of per^^A. 

Note that, as for determinant (see Section 2), for the receiving the re- 
quired enumerative information, we should prohibit to use the identities of 
type 1 + U + ... + cu™-^ = 0. 

2) Cycl{A). For n > 1, here we should ignore element an. Consider all 
full cycles vr with the condition 7r{j) = 1. Let j corresponds to j-th column 
of matrix A. Then the considered full cycles correspond to diagonals of ma- 
trix A having the common position (1, j), j >'2. Consider continuation of 
these diagonals in the matrix of combinatorial miner Aij. Then, by Lemma 
[1], the number of cycles of every its diagonal does not change, i.e., they are 
full cycles of of order n — 1. Therefore, we have the following expansion of 
Cycle{A) over the first row of A 

n 

(4.2) Cyde{A) = Y,aiACM),„ 

where (CM)ij, j > 1, are combinatorial minors of Cycle{A). 

3) S{A;n,k). From very close to 1) arguments, we have the following 
expansion of S{A; n, k) over the first row of A 

n 

(4.3) S{A-n,k) = a,,,{CM)[\-'^ + ^ai,,(CM)g\ 

i=2 

where {CM)[f, J > 1, are combinatorial minors of S{A] n, k). 
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Note that close to (4.3) formula was found by the author in [6] but using 
much more complicated way. 

4) C{ti,t2, ...jtn). We need lemma. 
Lemma 2. Let 

(4.4) {ctij! c^fcii, Ofcafci, ak^,kr-ii o-iMr) 
he a cycle. Then 

(4.5) {flfcll, 0^fc2fcl, «fcr,fer-ll '^i.fcrl 

is a cycle with respect to the main diagonal of the matrix of combinatorial 
miner Aij. 

Proof. According to the construction of Aij, the main its diagonal is 

(4.6) {022, ^33) •••) %1) '^i+lj + l) •••) O^nn}- 

With respect to this diagonal we have the following contour which shows 
that fl4.5p is, indeed, a cycle. 

(4.7) {ofe-^i — )■ aji — )■ Ojk^ — )■ ak,.kr-i O'kr-ikr-2 ~^ ■■■ ~^ flfc2A.i(~^ cifcii)}- 
■ 

Quite analogously we can prove that to every another cycle of a diagonal 
containing element aij correspond the same cycle with respect to the main 
diagonal of the matrix Aij. 

Let A be (0,1) square matrix of order n. Denote C^'^\A]ti,t2, ...^tn) a 
partial cyclic indicator of indicator (11. 9p of permutations vr G B{A) for 
which {1, 7r(l), 7r(2), 7r(r — 1)} is a cycle of length r. Then we have 

n 

(4.8) Yl ^^'^(^5 ^1' ^2, tn) = C{A; t,, t2, t„). 

r=l 

Therefore, it is sufficient to give an expansion of C^^^{A]ti,t2, ...jtn), r = 
1, n. First of all, note that 

(4.9) C«(A ^1, t2, tn) = antiCiA; tiM, tn-i). 
Furthermore, using Lemmas [T]121 we have 

n 

(4.10) C^'-\A-MM,-,tn) = Vai,,(CM)i,„ 
where 

(4.11) (CM)i,, = C(^-^)(Ai,,;ti,t2,...,tn-i), J = 2,...,ri, 

are the combinatorial miners of cyclic indicator of permutations with re- 
stricted positions. 
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Note that factor in ( l4.1Up corresponds to the diminution of the length 
of cycle f l4.5p with respect to length of cycle (14. 4p . Thus formulas (14. 8p - 
(4.11) reduce calculation of cyclic indicator of n-permutations to n — 1- 
permutations with a rather simple computer realization of this procedure. 
Note that similar but much more complicated procedure was indicated by 
the author in [7]. 

5. An example of enumerating the permutations of classes 
0,1,2 modulo 3 with restricted positions 



Let 



A 



(\ 

1 
1 



o\ 
1 

1 



V 



Consider class Bi^A) of permutations with restricted positions and find the 
distribution of them over classes 0,1,2 modulo 3. We use w-permanent 
with u) = e~ and its expansion over elements of the first row, given by 
(14. ip . Recall that, for the receiving the required enumerative information, 
we cannot use identities of type 1 + u + u'^ = 0. 
We have 

/l 1 l\ 

111 



per, 




oj^per 




{u + u^) + {uj + u'^ + l) + {2u^ + 1) + 



[uj- 



+ 2 + w) + (2 + 2a;) + (w^ + I) + {u'^ + I + u) + {2 + u) + 
{u + 2a;2 + 1) + (a;2 + 2 + w) = 13 + 9w + lOw^. 
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Thus in B{A) we have 13 permutation of class 0; 9 permutations of class 
1 and 10 ones from class 2 modulo 3. 



6. On TWO SEQUENCES CONNECTED WITH Cycle{A) 

In summer of 2010, the author published two sequences A179926 and 
A180026 in OEIS [8]. a{n) := Al79926(n) is defined as the number of per- 
mutations of the divisors of n of the form di = n, ^2, d^, dr{n) such 
that is a prime or 1/prime for all i. Note that a{n) is a function of 
exponents of prime power factorization of n only; moreover, it is invariant 
with respect to permutations of them. This sequence is equivalently defined 
as the number of ways, for a given finite multiset E, to get, beginning with 
E, all submultisets oi E, if in every step we remove or join 1 element of E. 
Sequence A180026 differs from A179926 by an additional condition: is 
a prime. In the equivalent formulation it corresponds to the condition that 
E is obtained from a submultiset in the last step, by joining 1 element. 

Note that, it is easy to prove that, knowing any permissible permutation 
of divisors, say, 5i = n, S2, Sr{n) (such that is a prime or 1/prime), 
we can calculate b{n) := 74l80026(n), using the following construction. Con- 
sider square (0,1) matrix B = {bij} of order T{n) in which bij = 1, if |^ is 
prime or 1/prime, and bij = 0, otherwise. Then b{n) = Cycle{B). In case of 
A179926, the construction is a little more complicated: a{n) = Cycle{A), 
where A is obtained from B by the replacing the first its column by the 
column from I's. 

Note also that A. Heinz [2] proved that, in particular, a{Il2^iPi), where 
Pi are distinct primes, equals to the number of Hamiltonian paths (or Gray 
codes) on n-cube with a marked starting node (see A003043 in [8]), while 
b{Ilf^iPi) equals to the number of directed Hamiltonian cycles on n-cube 
(see A003042 in [8j). 
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